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We compute the QCD phase diagram in the plane of the chiral chemical potential and temperature 
using the linear sigma model coupled to quarks and to the Polyakov loop. The chiral chemical 
potential accounts for effects of imbalanced chirality due to QCD sphaleron transitions which may 
emerge in heavy-ion collisions. We found three effects caused by the chiral chemical potential: 
the imbalanced chirality (i) tightens the link between deconfinement and chiral phase transitions; 
(ii) lowers the common critical temperature; (iii) strengthens the order of the phase transition by 
converting the crossover into the strong first order phase transition passing via the second order 
end-point. Since the fermionic determinant with the chiral chemical potential has no sign problem, 
the chirally imbalanced QCD matter can be studied in numerical lattice simulations. 
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I. INTRODUCTION 

It is well known that the QCD vacuum has a non- 
trivial topological structure due to presence of certain 
ghion configurations, instantons, which are characterized 
by an integer- valued topological winding number p!] . For 
a long time an experimental evidence for the existence of 
the topological gluon configurations could only be found 
indirectly, in certain features of the meson spectrum [2]. 

Recently, it was noticed that a potentially observable 
direct signature of the topologically nontrivial gluon con- 
figurations can emerge in noncentral collisions of heavy 
ions [31 14] . Such collisions create hot expanding fireballs 
of the quark-gluon plasma in the background of a strong 
magnetic field. Topologically nontrivial sphaleron transi- 
tions E] can induce - acting via the axial anomaly - a 
nonvanishing chiral density of quarks in the plasma fire- 
balls. This chirally imbalanced matter is characterized 
by different densities of right- and left-handed quarks. If 
such media is placed in an external magnetic field, then 
an electric current of quarks should emerge along the 
field's axis. The unusual generation of the electric cur- 
rent in the magnetic field background is the essence of 
the so-called "chiral magnetic effect" [3J H] . Signatures 
of this phenomenon were found in heavy-ion experiments 
at the BNL Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven [7] and they may also be observed in the 
heavy-ion collisions at the Large Hadron Collider (LHC) 
at CERN. 

The chiral magnetic effect is realized in the chirally im- 
balanced media in the strong magnetic field background. 
It is also important that the system is hot enough to 
be in the deconfining and chirally restored phase. The 
later requirement is essential because the system should 
be able to generate the quark's electric current so that 
the quarks should be deconfined. Moreover, in the chi- 
rally broken (low temperature) phase the chiral imbal- 



ance should quickly be washed out due to the presence 
of the chiral condensate which facilitates transitions be- 
tween left and right quarks. Therefore, it becomes inter- 
esting to investigate the infiuence of the background mag- 
netic field and the effect of the chiral imbalance on the 
finite-temperature transition between hadron and quark- 
giuon plasma phases. 

The effect of the strong magnetic field on the QCD 
phase transition was studied both analytically [8l-fTT] and 
numerically [T^. All these studies have found that the 
magnetic field background increases the transition tem- 
perature and makes the phase transition stronger. A 
third effect was found in Ref. [10]: the magnetic field 
splits the deconfining and chiral phase transitions, thus 
leading to emergence of the new, chirally broken decon- 
fining phase (the splitting can be small [TT], however). In 
the low-temperature and strong-magnetic-field corner of 
the QCD phase diagram a new electromagetically super- 
conducting phase may emerge |13'. 

The topologically induced changes in chirality can be 
modeled with the help of the chiral chemical potential 
/is which creates a difference between the right- and left- 
handed particles [31 . This potential can be related to the 
6 angle of strong interactions as follows [4] : 

" dt 2Nf ' ^ ' 

where t is the time coordinate and Nf is the number of 
the light flavors in the theory. 

The influence of the chirally imbalance on the ther- 
mal phase transition in the magnetic field background 
was addressed in Ref. [9^ regarding possible applications 
to the chiral magnetic effect. Working in the Nambu- 
Jona-Lasinio (NJL) model coupled to the Polyakov loop 
(PNJL), the authors of Ref. J] have found that the chi- 
ral imbalance makes the temperature of the chiral phase 
transition smaller while the strength of the transition 
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becomes stronger. In our paper we confirm findings 
of Ref. [9] working at zero magnetic field in the linear 
sigma model coupled to quarks and to the Polyakov loop 
(PLSMq) which also serves as an effective low-energy 
model of QCD. Our main result is the QCD phase di- 
agram in the (/^5, T) plane, which is plotted in Fig. |4] 

The structure of this paper is as follows: in Section |TT] 
we describe our model (PLSMq), in Section III we discuss 
thermodynamics and calculate the phase diagram of the 
model, and Section [IV| is devoted to our conclusions. 



II. THE MODEL 

We use a linear sigma model coupled to quarks [Ti] 
and to the Polyakov loop (PLSM^) [15]. The inclusion of 
the Polyakov loop allows us to account for effect of the 
color confinement following similar proposal in the NJL 
model [16j. 

This low-energy model of QCD contains three types of 
fields: the duplet of the quark fields ip{x) = {u,d)^, the 
scalar chiral (meson) fields (ct, 7?) with 7? — (tti, 7r2, tts), 
and the complex-valued scalar field of the Polyakov loop. 



L{x) = -Tr^x) 



$ = P exp 



l/T 



dTA4x,T) ,(2) 



where A4 — iAo is the timclike component of the SU (3) 
gauge field, V is the path-ordering operator, if is the 
isotriplet of the pseudoscalar pions, and a is the pseu- 
doscalar field. 

The PLSMg Lagrangian can be represented as a sum 
of the following three parts: 



C = Cq{ip,ip,(T,Tf,L) + £o-(ct, tt) + Cl{L) 
The quark part of the Lagrangian ([s]). 



ijj [ilp - g{a + i^f ■ tt) -f MsT 7 V' 



(3) 



(4) 



provides the interaction between the quarks %p, the chi- 
ral fields (J, TT, and the gauge field via the covariant 
derivative Ip — 7^(9^ — iA^). The Lagrangian Q in- 
cludes also the real- valued chiral chemical potential /is. 

The dynamics of the chiral fields is described by the 
second term in the Lagrangian (|3]), 

-K(^,Tr), (5) 
where we have introduced charged and neutral mesons. 
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^0 ?, 

vr = vr . 



(6) 



respectively. The first and the second terms of the po- 
tential, 



v^) — ha , 



(7) 



provide, respectively, (strong) spontaneous and (weak) 
explicit breaking of the chiral symmetry. The phe- 
nomenologically acceptable parameters in Eqs. Q and 
^ are: g = 3.3, A = 20 and w = 87.7 MeV PT!. We work 
in a mean field approximation thus neglecting quantum 
fluctuations (i.e., the kinetic terms) of the scalar fields a 
and TT. 

Effects of the color conflnement are encoded in the last 
term of the Lagrangian ([3| , which describes the potential 
of the Polyakov loop: 



-VUL, T) 



(8) 



As in the previous case involving the mesonic chiral fields, 
we neglect possible kinetic terms for the Polyakov loop 
and consider only the potential term following Ref. |18) : 



Vl{L,T) ^ 1 
T4 2 



a{T)L*L 



(9) 



fe(r)ln 1 -6L*L + a(l*^ + L^) -3{L*Lf , 



where 



a{T) = Oo + a 
'To 



h{T) = 63 



To 
1 I ^ 

3 



a2 



T 



(10) 



and To = 270 MeV is the temperature of the deconfine- 
ment phase transition in the pure SU{i) Yang-Mills the- 
ory without quarks. The coefficients in Eq. ( 10 1 are: 



ao = 16 7r745 « 3.51, ai 
a2 = 15.2, 63 



-2.47, 
-1.75. 



(11) 



The Polyakov loop potential (|9| respects the center Z3 
symmetry, L — >• e^'^"'/^L with n = 0, 1, 2. 



III. THERMODYNAMICS AND PHASE 
DIAGRAM 

The free energy density of the system is: 



T 

17 = -— InZ, 



(12) 



where Z is the partition function and V^d is the volume 
of the three-dimensional space. In the mean field ap- 
proximation the chiral fields and Polyakov loop field are 
considered as classical objects, so that the free energy 
can be rewritten as 

n (a, 7f, L, T, /X5) = V^{a, tt) + Vl {L, T) 

+^q{a, 7?, L, T, /ig) , (13) 

where the potential is given by Eq. ([t]) and the poten- 
tial Vl is described in Eqs. (|9|-([Tl]). In Eq. ([l3]) the quark 
part flq comes from the fermion determinant, which can 
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be computed by explicit diagonalization of the quadratic 
operator corresponding to the quark Lagrangian Q: 



-2E 



=±1 



(27r) 



(14) 



-r|in[i 



1+3 iL*+Le' 



s/T 



+C.C. 



where L is the Polyakov loop ([2| and the fermion spec- 
trum is 



ujsip) = v/(b|s-M5)2+5'(f^'+^') 



(15) 



Here s = ±1 is the helicity (the sign of the projection 
of the particle's spin on the direction of the particle's 
motion). The first term in Eq. (14 1 corresponds to the 



divergent energy of the Dirac sea. After a proper regu- 
larization the contribution from this term renormalizes 
the parameters A and v of the chiral potential ([7| adding 
also finite logarithmic corrections (which are sometimes 
called zero-point corrections) to the potential. These cor- 
rections have a pure vacuum origin because they depend 
neither on temperature nor on the chiral chemical po- 
tential. Usually, these logarithmic corrections are disre- 
garded in the phenomenological approaches based on the 
linear sigma model because their contribution does not 
change the qualitative physical picture. We ignore the 
logarithmic corrections as well, referring the interested 
reader to studies of the vacuum fluctuations effects done 
in Refs. [T1I2D]. 



Integrating Eq. ( 14 ) by parts, we can rewrite the free 



energy in the following form: 



s=±l 



..duJsip) , , 

^9,^(P)J^^'(16) 



where nq^s and n^.s = n* ^ are, respectively, the occupa- 
tion numbers (summed over colors) for quarks and anti- 
quarks carrying the spirality s, and 



3e-"=/^ (L + 2L*e-"=/^ + g-^-./T') 
1 + 3 (L + L*e--''/T) e-^'/T + g-s^,/! 



(17) 



The mean field values of the fields ct, t? and L are 
found by a (numerical) minimization of the free energy 
given by Eqs. ([l3|, Q, (|9|-([lT]), and ([l6]) with respect 
to the variations of these fields at fixed values of the 
temperature T and chiral chemical potential ^l^. 

The fermion determinant breaks the center Z3 sym- 
metry of the Polyakov loop potential and the global min- 
imum of the free energy ( 13 ) corresponds to the real- 
valued Polyakov loop. Moreover, one can show numeri- 
cally that the presence of the charged pion condensates 
TT^ makes the free energy larger so that these conden- 
sates are disfavored, {-k^) — 0. The last term in Eq. Q 
forces the neutral pion condensate to be zero, (7r°) — 0. 

Therefore we are left with two unknown expectation 
values which are to be fixed by the minimization of the 
free energy. These are the real part of the Polyakov loop, 
ReL, and the chiral order parameter <t. 



The expectation values of the Polyakov loop and the 
chiral order parameter a are shown as functions of the 
temperature in Fig. [T] and Fig. [2] respectively. The chiral 
field is normalized to unity at T = /is = 0, and the 
corresponding vacuum expectation value ctq is given by 
the pion decay constant: 



CTO 



(a)T= 



f^{T^^i^^O) = 92.2 MeV . (18) 



The temperature and chiral chemical potential in Figs. [T] 
and [2] are expressed in units of the critical transition tem- 
perature Ti"-* at zero chiral chemical potential in PLSM^: 

T^(°) EE Tc(/i5 = 0) = 213.2 MeV [in PLSM,] . (19) 

Due to the presence of the explicit symmetry breaking 
term in the chiral potential Eq. ([t]), the transition 
between the quark-gluon plasma phase and the hadron 
phase at zero chemical potential is a smooth crossover. 
The critical temperature value in Eq. ( jig] ) corresponds to 
the temperature where the slopes of the Polyakov loop L 
and the chiral field a are steepeslj^ 
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FIG. 1. The expectation value of the Polyakov loop as the 
function of temperature T at fixed values of the chiral chem- 
ical potential ^5 [T and /is are given in units of ri"\ the 



critical temperature of the fis — transition, Eq. ( 19 1 



According to Figs.[T]and[2j the crossover turns into the 
first order transition as the chiral chemical potential /is 
increases (the larger /is , the stronger the transition) . The 
critical temperature of the transition between the quark- 
gluon plasma phase and the hadron phase decreases as 
the chiral chemical potential increases. Moreover, the 
chiral and deconfinement phase transitions are tightened 
to each other and they do not split as the value of /is 
increases. 



^ The difference in the critical values for these two quantities is 
smaller then one MeV in PLSM„. 
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FIG. 2. The chiral order parameter a (we use the same nota- 
tions as in Fis. The value of £70 is given in Eq. ||18[). 



Obviously, the presence of the nonzero chiral chemi- 
cal potential /15 does not induce a nonvanishing baryon 
charge, so that the quark density n is always zero. 



s=±l 



{2ny 



[nqAp) 



'■g,s(p)J 



0. (20) 



However, the density of the chiral charge (i.e., the differ- 
ence between the densities of the right- and left-handed 
particles). 



n5 =nR- riL 



3d 



(21) 



should in general be nonzero at a nonvanishing chiral 
chemical potential (and, naturally, 77.5 = at /is = 0). 

We plot (the one-third power of) the density of the 
chiral charge ( [2I] ) is Fig. [s] The chiral imbalance of the 
media in the presence of a fixed chemical potential /is 
is a growing function of the temperature. The chiral 
charge density ns is strongly enhanced at the transition 
temperature Tc — Tdp.^), as the system goes from the 
hadron phase to the quark-gluon plasma phase. 

The phase diagram in the (/is , T) plane is shown in 
Fig. [4] The transformation of the crossover transition 
(the dashed line at lower values of /is) to the first order 
phase transition (the solid line at larger /is) goes via a 
critical endpoint (CEP) in which the transition becomes 
a second order transition. We found that the CEP (the 
red point in Fig. ^ is located at 

(/ts,T)^^^ = (0.232 ri"), 0.998 Ti°)) . (22) 

Our phase diagram, Fig. |4j computed in the PLSM, 
at zero magnetic field agrees qualitatively with the cor- 
responding phase diagram of Ref . [9^ , where the calcula- 
tions are done in the PNJL model in a weak magnetic 
field backgroun(|^ in the window < /^s < 2.3ri"\ The 



^ The critical temperature curve Tc = Tc(/i5, B) is almost indepen- 
dent on the strength of the magnetic field B in a weak field [9]. 



FIG. 3. The one-third power of the chiral charge density n^, 
Eq. as a function of temperature T at a fixed set of 

the chiral chemical potential /15 (we use the same units as in 

Fig.g. 
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FIG. 4. Phase diagram of QCD in the plane of the chiral 
chemical potential /is and temperature T [expressed in uni ts 



of the critical transition temperature ri''' at /is 
The dashed line represents the smooth crossover 
line corresponds to the first order phase transition, and the 
red point marks the second-order critical endpoint (22 I. The 



0,Eq. figf]. 
the solid 



inset shows the same phase diagram with the chiral chemical 
potential axis plotted in the logarithmic scale. 



critical curves = Tc{fi5) in the PLSM^ and PNJL mod- 
els are quite close to each other, while the positions of 
the corresponding CEP's differ from each other substan- 
tially, as the location of the CEP in the PNJL model 



is (/i5,r) 



CEP _ 



(1.7ri°\o.8Ti''0 according to Ref. [9]. 



The difference between the PNJL prediction and our re- 



sult ( 22 ) can be attributed either to the mo del- dependent 
issues or to the fact that we have neglected the vac- 
uum corrections coming from the fermionic determinant 
which may influence the PLSMg thermodynamics quan- 
titatively [20]. 

It is important to mention that the line of the first 
order phase transition does not hit the /is axis at any 
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finite value of the cliiral cfiemical potential. Instead, the 
PLSMg model predicts that the Tc = Tdii^) curve ap- 
proaches this axis smoothly at asymptotically large val- 
ues of the chiral chemical potential /is : 



lim Tcifis) = . 



(23) 



Notice, however, that PLSM, is an effective low-energy 
model of QCD, and therefore at large values of massive 
parameters (for example, at large fj,^) the results coming 
from this model may become inaccurate. Therefore the 



prediction in Eq. (23) should be considered with care. 



IV. CONCLUSIONS 

The chirally imbalanccd hot quark-gluon plasma may 
emerge in heavy-ion collisions at RHIC and LHC exper- 
imental facilities. We have computed the QCD phase 
diagram in the chirally imbalanced background at finite 
temperature (Fig. [4| using the linear sigma model cou- 
pled to quarks and to the Polyakov loop (FLSM^) at zero 
magnetic field. Our results are in a qualitative agree- 
ment with conclusions of an earlier study of the QCD 
phase finite-temperature phase transition in the PNJL 
model in simultaneously imposed chiral and magnetic- 
field backgrounds [9]. 

We have found that the increase of the chiral chemi- 
cal potential /is tightens the link between deconfinement 
and chiral phase transitions, simultaneously lowering the 
common critical temperature and strengthening the or- 
der of the phase transition by converting the crossover 
(realized at low chiral imbalance) into the first order 
phase transition (found at higher chiral imbalance) . The 
location of the second-order critical endpoint - at which 
the crossover turns into the first order phase transition - 



is given in Eq. ( 22 ) . The PLSMg predicts that at strictly 
zero temperature the system stays always in the hadronic 
phase regardless of the value of the chiral chemical po- 
tential [Eq. ([23]) and Fig.|4]. 

We are mainly interested in the finite-temperature 
phase diagram of QCD at nonzero chiral chemical po- 
tential because the chiral magnetic effect is realized in 
the presence of the strong magnetic field in the chirally 
imbalanced background [Ij. We confirm that the some 
effects of the strong magnetic field and the chiral chem- 
ical potential on the critical temperature of the QCD 
phase transition are opposite: the magnetic field tends 
to increase the temperature of the phase transition [S]- 
[T2| while the chirally imbalanced background forces the 
critical temperature to become lower 9J. Moreover, the 
magnetic field background splits the chiral and deconfine- 
ment phase transitions [10] (the splitting can be small, 
however [H]), while the chiral imbalance tightens the 
link between these transitions. However, despite these 
dissimilarities, the external magnetic field and the chiral 
imbalance have one common feature: they both make the 
QCD phase transition stronger. 

Finally, we would like to notice that the predictions 
for the QCD phase diagram in the (/xs, T) plane, given in 
Fig.|4]of the present article and in Ref. [9], can directly be 
tested from the first principles in numerical simulations 
of lattice QCD because the fermionic determinant with 
the chiral chemical potential has no sign problem [3]. 
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